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Abstract. Path integral techniques in collective fields are shown to be a useful ana- 
lytical tool to reformulate a field theory defined in terms of microscopic quark (gluon) 
degrees of freedom as an effective theory of collective boson (meson) fields. For illus- 
trations, the path integral bosonization approach is applied to derive a (non)linear a 
model from a Nambu-Jona-Lasinio (NJL) quark model. The method can be extended 
to include higher order derivative terms in meson fields or heavy-quark symmetries. It 
is also approximately applicable to QCD. 

1 Introduction 

In this lecture, I want to demonstrate the powerfulness of the path integral 
approach in collective fields for the bosonization of quark models containing 
(effective) 4-quark interactions [1]. For illustrations, let me consider NJL type of 
models [2-4] with local quark interactions ~ G {q_rq) representing a relativistic 
version of the superconductor BCS theory [5,6]. These models lead to a gap 
equation for a dynamical quark mass signalling the dynamical breakdown of 
chiral symmetry (DBCS). Furthermore, the collective field of Cooper pairs of 
the superconductor is now replaced by collective meson fields of (gg)-bound 
states. 

To be more explicit, let me consider the generating functional Z of the NJL 
model defined by a path integral of the exponential of the corresponding action 
over quark fields q as the underlying microscopic degrees of freedom. Path in- 
tegral bosonization then means to transform this generating functional into an 
integral of the exponential of an effective meson action where the new collective 
integration variables a, 7r, p, ... denote the observable meson fields, 

Z = J Dfi(q)e [ I CN1Ld4x =L> J D t i(a,ir,p,...)e i S Cmt - A * x , (1) 

with D/i(q) = DqDq, Dfx (a, 7r, p, ...) = DaDirDp... being the respective inte- 
gration measures of fields. 
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The basic ingredient of the path integral bosonization (1) is the use of the 
Hubbard-Stratonovich transformation [7,8] which replaces the (effective) 4-quark 
interactions of NJL models by a Yukawa-type coupling of quarks with collective 
meson fields fa — (<r, 7r, p„, ...). After this the primary path integral over quark 
fields on the L.H.S. of (1) becomes Gaussian resulting in a quark determinant 
containing meson fields. Further important steps are: 

(a) the use of the loop expansion of the quark determinant in powers of 
meson fields 

(b) the evaluation of the resulting Feynman diagrams in the low-momentum 
region 

(c) the limit of large numbers of colours, N c — > oo, GN C fixed, G being the 4- 
quark coupling constant in order to apply a saddle point approximation 
to the integration over meson fields. 

Quark loop diagrams emitting two ^-fields then generate in the low-momentum 
(two-derivative) approximation kinetic and mass terms of mesons. Finally, quark 
loops emitting n > 2 meson fields lead to meson interactions with effective small 

coupling constants g n <~ O ^(y^r) ^ allowing for a modified perturbation 

theory in terms of meson degrees of freedom. 

Notice that the NJL quark Lagrangian incorporates the global chiral flavour 
symmetry of Quantum Chromodynamics (QCD) as well as its explicit and dy- 
namical breaking. The equivalent effective meson theory on the R.H.S. of (1) just 
reproduces this symmetry breaking pattern at the meson level. As mentioned 
above, masses and coupling constants of collective mesons are now calculable 
from quark loop diagrams and expressed by the parameters of the underlying 
quark theory (including a loop momentum cut-off A). The final aim is, of course, 
bosonization of QCD, i.e. the transformation 



Z = J DufaGjjf**™ 6 '*^ J D^{a^,p,...)e l I c ^ x , (2) 

where q, G M denote quark and gluon fields (ghost fields are not shown explicitly). 
In order to begin with an effective 4-quark interaction as intermediate step, one 
first has to integrate away the gluon (ghost) fields on the L.H.S. of (2). This can 
only be done exactly for space-time dimensions D = 2 in the light-cone gauge, 
G- = \ (Go — G\) = 0, where all self-interactions of gluon fields vanish. The 
resulting expression contains an effective nonlocal current x current quark interac- 
tion with a known gluon propagator which then can be bosonized by introducing 
bilocal meson fields (j) (x, y) ~ q(x)rq(y) 1 employing the limit of large numbers of 
colours iV c [9]. Clearly, the analogous bosonization of four-dimensional QCD is 
more complicate. It requires some additional approximations: first a truncation 
of multilocal quark current interactions retaining only the bilocal two-current 
term and secondly a modelling of the unknown nonperturbative gluon propaga- 
tor (Cf. Fig. la) (see [1,10-13]). As a final result one derives an effective chiral 
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Lagrangian describing the low-energy interactions of observable mesons includ- 
ing nontrivial meson form factors. 




a) b) 

Fig. 1. a-b) Low-energy approximation of a nonlocal current x current interaction with 
nonperturbative gluon propagator (a), by a local Nambu-Jona-Lasinio type interaction 
(b). 

Nevertheless, as anticipated in (1), also simpler local NJL type of models 
are expected to yield a reasonable low-energy description of the chiral sector 
of QCD 4 . Indeed, discarding the complicated question of the unknown struc- 
ture of nonperturbative gluon and quark propagators (related to confinement), 
one can try to approximate the nonperturbative gluon propagator in the re- 
gion of low energies by a universal constant G leading to a local NJL type of 
current x current interaction (Cf. Fig. 1). (For interesting nonperturbative ap- 
plications of the above approach the reader is further referred to the nuclear 
many-body problem [10,14] and the Hubbard model [15,16].) 

In conclusion, let us notice that for two space-time dimensions there exists 
a different realization of the bosonization idea based on an explicit construction 
of fermionic fields in terms of bosonic fields due to Mandelstam [17]. This then 
allows one to replace, at the operator level, quark bilinears by bosonic fields, e.g. 

ipltiip ~ -n~^£iivd v (p, tpip ~ M cos (2y/n(j)) , 

and to prove the equivalence of a given fermionic model (e.g. Thirring model) 
with a corresponding bosonic model (e.g. Sine-Gordon model). Using in partic- 
ular Witten's non-Abelian bosonization rules [18] one has derived in the strong 
coupling limit (which is contrary to the weak coupling limit, GN C fixed, for 
N c — > oo, of the above described path integral approach) a low-energy effective 
action from QCD 2 [19]. The lecture of T.Giamarchi at this Workshop [20] dis- 
cusses just the application of this kind of operator bosonization in Condensed 
Matter Physics. 
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In the next section I shall now describe the path integral bosonization of the 
NJL model along line / closely following the original papers [3,4]. 

2 NJL Model and cr Model 
2.1 Linear cr Model 

Let us consider the following NJL Lagrangian with a global symmetry 



where q denotes a quark spinor with colour and spinor indices, r,'s are the genera- 
tors of the flavour group U(2), and G is a universal quark coupling constant with 
dimension mass -2 . Note that the gg-combinations in (3) are colour singlets, and 
we have admitted an explicit chiral symmetry-breaking current quark mass term 
—m^qq. The integration over the quark fields in the generating functional Z of 
the NJL model (given by the L.H.S. of (1)) can easily be done after bi-linearizing 
the four-quark interaction terms with the help of colour-singlet collective meson 
fields a, 7r. To this end, we use the Hubbard-Stratonovich transformation [7,8] 



which replaces the 4-quark interaction by a Yukawa coupling with collective 
fields cr,7r, Cf. Fig.2. 



[U{2) x U(2)} x SU(N C ) [3] 




(3) 






Fig. 2. Graphical representation of the integral identity (4). 



Inserting (4) into the L.H.S. of (1) leads to the intermediate result 
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with the semi-bosonized meson-quark Lagrangian 

£ njl = - (O - m of +^ +q(id-cr- i 75 T • tt) q, (6) 

where we have absorbed the current quark mass in the a field by a shift a — > cr— m 
of the integration variable. The Gaussian integral over the Grassman variable q 
in (5) can easily be performed leading to the fermion determinant 

detS" 1 = expiV c Tr lnS^ 1 = cxpiiV c J -itr (in S" 1 ) {x x) d 4 x (7) 

where 



S 1 (x,y)= id x - a(x) - i-y 5 T ■ ir(x) 5 i (x-y) 



(8) 



is the inverse quark propagator in the presence of collective fields a, it and the 
trace tr in (7) runs over Dirac and isospin indices. Note that the factor A^ c results 
from the colour trace. Combining (5)-(7) we obtain 



Z = J DaDme' I CE!! ^' v)di 



(9) 



£es. (o-,7r) = -— ^ ((a - m f + n 2 ) -iN c tr In (id - a - vy 5 r ■ tt) . (10) 

Z(_r V /V / (x,x) 

Let us analyze (9) in the limit of large N c with GN C fixed where one can apply 
the saddle point approximation. The stationary point a = (a) = m, tv = 
satisfies the condition 



da 



0, (11) 



which takes the form of the well-known Hartree-Fock gap equation determining 
the dynamical quark mass m 

m = m + 8mGN c Ii(m), (12) 

A 

= (13 » 

with A being a momentum cut-off which has to be held finite (see Fig. 3). 

Note that in the chiral limit, mo — » 0, (12) admits two solutions m = or 
m 7^ in dependence on GN^ (GN c ) cHt . Thus, for large enough couplings we 
find a nonvanishing quark condensate (qq) ~ trS(x,x) signalling spontaneous 
breakdown of chiral symmetry. 

It is further convenient to perform a shift in the integration variables 
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m 



+ 




Fig. 3. Graphical representation of the gap equation (12). 



a — > to + 



7T 



In order to obtain from the nonlocal expression (10) a local effective meson 
Lagrangian we have to apply the following recipes: 

(a) loop expansion of the determinant in the fields a', 7r, i.e. 



N c tr In 



. | (id — mj 



1 



id -to VV^V, 



175 T- 



7T 

Wc 




+ 



finite diagrams. 
Fig. 4. Loop expansion of the fermion determinant. 



Here we have omitted the unimportant constant term N c tr In [id — mj 

and used the formula In (1 — x) = — x + ^- H . (The tadpole dia- 
gram linear in a' cancels by a corresponding linear term arising from the 
first term in (10) due to the gap equation (12).) 
(b) Low-momentum expansion of loop diagrams corresponding to the so- 
called gradient expansion of meson fields in configuration space. 
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(c) Field renormalization, 

(a',n)^zi (a',7r). 

Discarding the contributions of finite n-point diagrams with n > 4 in the loop 
expansion leads to a linear a model of composite fields (the prime in a' is now 
omitted) [3] 

£eaf. (ct, it) = iff (-□ -ml)a+ ^tv (-□ - m 2 ) tt- 

- 9o^O + TT 2 ) - 947T (CT 2 + TT 2 ) 2 . (14) 

Notice that the masses and coupling constants of mesons are fixed by the quark 
model parameters and the (finite) momentum cut-off A, 

m 2 =m ^g, to 2 =to 2 +4to 2 , (15) 



g* qq = g* qq = (JQ 2 - {4iv c / 2 } ' = o 



/Nr. 



8N C I 2 ■ 
where 



'o-tttt — 93a — T, (16) 

1 

.947T = <74<T 



h = -i 



d 4 k 



(2tt) 4 (fc 2 - to 2 ) 2 ' 



Notice that in the chiral limit too — ► 0, the pion becomes the massless Goldstone 
boson associated to DBCS. Introducing electroweak interactions at the quark 
level leads to the additional diagram shown in Fig. 5. 
Fig. 5 leads to the interaction Lagrangian 

A£= 9 -W;(-—d»A, (17) 

which yields just the PCAC meson current. Here the ratio m/g^qg has the mean- 
ing of the pion decay constant F n . We thus obtain the Goldberger-Treiman re- 
lation 

F« = —, (18) 

9ir qq 

valid at the quark level. 

In the following subsection we shall show how one can derive the related 
nonlinear version of the er model. 
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Fig. 5. Quark diagram describing the weak transition tt — ► W. 



2.2 Nonlinear er model 

Nonlinear chiral meson Lagrangians have been introduced a long time ago in the 
literature [21-23] and are now widely used in low-energy hadron physics [1,24]. 
Let us demonstrate how they can be derived from a NJL quark model. To this 
end, it is convenient to use in (6) an exponential parametrization of the meson 
fields 

a + i75T-7T = (m + a)e' h5 ^ , (19) 

which yields 

cr 2 + TT 2 = (m + af . (20) 

Let us absorb the exponential in (19) by introducing chirally rotated "con- 
stituent" quarks x, 

q = j'» z *x; £ = (21) 



leading to the Lagrangian 1 



+ X 



with 



16G 

being a symmetry-breaking mass term. 



i7 M (V + e" i75 a^ e i75 ^- m -a~] X , (22) 



ACsb = !! ^rn tr ( c -^< + h.c.) (23) 



1 There arises also an additional anomalous Wess-Zumino term from the fermion mea- 
sure of the path integral [4] which is discarded here. 
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Note the nonlinear transformation law of the meson field £(x) under global 
chiral transformations g G SU(2) A x SU(2) V [22,23], 

ge H5 z$P = e i75 h{x) , (24) 

where 

h(x)=jt- u M x M G SU(2) v ,ioc. 

is now an element of a local vector group. It is convenient to introduce the Cartan 
decomposition 

e-^V^ = i 75 ^ • A„ (€) + • V M (£) . (25) 

It is then easy to see that the fields transform under the local group SU(2)v,ioc. 
as follows (Afj, = § • Afj,, etc.) [23] 

x -» x = K x )x 

V M -» V; = h^V^x) - /i^i^/i^a;) 

^^^/i^-V^W- (26) 

Thus, is a gauge field with respect to SU(2)v,ioc- This allows one to define 
the following chiral-covariant derivative of the quark field \ 

D„x=(d fl +iV fl )x- (27) 

Using (22), (25) and (27), the inverse propagator of the rotated \ ne ld takes the 
form 

S^ 1 = iD-m-a-Aj 5 . (28) 

The nonlinear a model is now obtained by "freezing" the a field, performing the 
path integral over the x ne ld and then using again the loop expansion for the 
quark determinant detS^ 1 . Choosing a gauge- invariant regularization, the loop 
diagram with two external A^ fields generates a mass-like term for the axial 
Af,. (£) field contributing to the effective Lagrangian 2 

^ lin . = ^tr F Al + AC sb . (29) 

9-jrqq 

Here, the symmetry-breaking term is taken over from (22), and the first constant 
term has been omitted. 

2 Note that usually generated field strength terms —\A^ v ,—\Vf lv vanish for Car- 
tan fields (25). A mass-like term ~ V 2 does not appear due to gauge-invariant 
regularization. 
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As has been shown in [23], A.^ is just the chiral-covariant derivative of the £ 
field admitting the expansion 

4 = D,e = d,e + o (e) = ^ v + o . m 

Thus, we obtain the nonlinear a model 

-Cniin. - tfDpCD"? + AC sh = (31) 

= lv{-n-ml)<p + 0(<f>l) 

reproducing Weinberg's result [21]. 
3 Conclusions and Outlook 

In this talk I have shown you that the path integral bosonization approach 
applied to QCD-motivated NJL models is a powerful tool in order to derive 
low-energy effective meson Lagrangians corresponding to the (nonperturbative) 
chiral sector of QCD. 

The above considerations have further been extended to calculate higher- 
order derivative terms in meson fields by applying heat-kernel techniques to the 
evaluation of the quark determinant [4]. This allows, in particular, to estimate 
all the low-energy structure constants Li introduced by Gasser and Leutwyler 
[24]. Moreover, it is not difficult to include vector and axial- vector currents into 
the NJL model and to consider the chiral group U(3) x U(3). In a series of 
papers [1,3,4,25] it has been shown that the low-energy dynamics of light pseu- 
doscalar, vector and axial-vector mesons is described surprisingly well by effec- 
tive chiral Lagrangians resulting from the bosonization of QCD-motivated NJL 
models. These Lagrangians embody the soft-pion theorems, vector dominance, 
Goldberger-Treiman and KSFR relations and the integrated chiral anomaly. 

Finally, we have investigated the path integral bosonization of an extended 
NJL model including DBCS of light quarks and heavy quark symmetries of heavy 
quarks [26] (see also [27]). This enables one to derive an effective Lagrangian of 
pseudoscalar, vector and axial-vector D or B mesons interacting with light ir, p 
and ai mesons. Note that the use of the low-momentum expansion in the eval- 
uation of the quark determinant restricts here the applicability of the resulting 
effective Lagrangian to such (decay) processes where the momentum of the light 
mesons is relatively small. 

Summarizing, I thus hope to have convinced you that the path integral 
bosonization approach to QCD (in its bilocal or simpler local formulation) is 
a very interesting nonperturbative analytical method which, being complemen- 
tary to numerical studies of lattice QCD, is worth to be developed further. In the 
next lecture, I will consider the related path integral "hadronization" of baryons. 
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